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ABSTRACT 

Motivated by the recent results of Hansen et al. (2008) concerning a noticeable hemispherical power 
asymmetry in the WMAP data on small angular scales, we revisit the dipole modulated signal model 
introduced by Gordon et al. (2005). This model assumes that the true CMB signal consists of a 
Gaussian isotropic random field modulated by a dipole, and is characterized by an overall modulation 
amplitude, A, and a preferred direction, p. Previous analyses of this model has been restricted to very 
low resolution (ie., 3.6° pixels, a smoothing scale of 9° FWHM and £ < 40) due to computational cost. 
In this paper, we double the angular resolution (ie., 1.8° pixels and 4.5° FWHM smoothing scale), 
and compute the full corresponding posterior distribution for the 5-year WMAP data. The results 
from our analysis are the following: The best-fit modulation amplitude for £ < 64 and the ILC data 
with the WMAP KQ85 sky cut is A = 0.072 ± 0.022, non-zero at 3.3cr, and the preferred direction 
points toward Galactic coordinates (I, b) = (224°, —22°) ± 24°. The corresponding results for £ < 40 
from earlier analyses was A = 0.11 ± 0.04 and (l,b) = (225°,— 27°). The statistical significance of 
a non-zero amplitude thus increases from 2.8c to 3.3c when increasing £ max from 40 to 64, and all 
results are consistent to within la. Similarly, the Bayesian log-evidence difference with respect to 
the isotropic model increases from A\nE — 1.8 to Ami? = 2.6, ranking as "strong evidence" on the 
Jeffreys' scale. The raw best-fit log-likelihood difference increases from Aln£ = 6.1 to Aln£ = 7.3. 
Similar, and often slightly stronger, results are found for other data combinations. Thus, we find that 
the evidence for a dipole power distribution in the WMAP data increases with £ in the 5-year WMAP 
data set, in agreement with the reports of Hansen et al. (2008). 

Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 



1. INTRODUCTION 

The question of statistical isotropy in the cosmic mi- 
crowave background (CMB) has received much attention 
within the cosmological community ever since the release 
of the first-year Wilkinson Microwave Anisotropy Probe 
(WMAP; Bennett et al. 2003a) in 2003. The reasons for 
this are two-fold. On the one hand, the current cosmo- 
logical concordance mode l is based on the concept of in- 
flation (lStarobinsky||1980l: iGuth et al II 19811: [Liride et al l 
19821: iMukhanov et al.l 119811: IStarobinskv et all 119821: 
Linde et all 119831. 119941: ISmoot et all 119921 : iR uhl et al I 
20031: iRvnvan et all 120031 : IScott et al.ll2003D . which pre- 
dicts a statistically homogeneous and isotropic universe. 
Since inflation has proved highly successful in describ- 
ing a host of cosmological probes, most importantly the 
CMB and large-scale power spectra, this undeniably im- 
poses a strong theoretical prior towards isotropy and ho- 
mogeneity. 

On the other hand, many detailed studies of the 
WMAP sky maps, employing higher-order statistics, 



have revealed strong hints of both violation of statis- 
tical isotropy and non-Gaussianity. Some early no- 
table examples include unexp ected low-l correlations 
(jde Oliveira-Costa et al.l l2004h . a peculia r large cold 
spot in the southern Galactic hemisphere (jVielva et al.l 
|2004|). and a dipolar distribution of large-scale power 
( Eriksen et al.l l2004af) . Today, the literature on non- 
Gaussianity and violation of statistical isotropy in 
the WMAP data h as grown very large, i ndeed (e.g., 
Bernui et alj 120061: iBielewicz et all 120051: ICopi et al.l 



2006UCruz et al.l200 5. 2006: Eriksen et al . 2004b. c EoolT 
Jaffe et al.1 120051 120061: iMartmez-Gonzalez etail 120061: 

Wandclt 



McEwen et al.ll2008l : iRath et al.ll2007HYadav _. 
20081) . and it would be unwise not to consider these issues 
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very seriously. 

Of particular interest to us is the question of hemi- 
spherical dis tribution of power in t he WMAP data, first 
reported by Eriks en~et al.l (l2004al) and later confirme d 
by, e.g., lHansen et all (|2004l ) and [Eriksen et al.l (|2005h . 
The mos t recent works on th is topic include those pre- 
sented bv lHansen et al.l (|2008f ) , who found that the power 
asymmetry extends t o much smalle r scales than previ- 
ously thought, and bv lEriksen et al.l (|2007bl ). who quan- 
tified the large-scale power asymmetry in the 3-year 
WMAP data using an optimal Bayesian framework. 

A separate, but possibly physicall y related, line of 
work was recently presented by iGroeneboom &: Eriksen! 
(2009), who considered the specific model for viola- 
tion of L orenz invariance in th e early universe, pro- 
posed bv lAckerman et ail (|2007h . This model involves 
CMB correlations with a quadrupolar distribution on 
the sky, and is thus orthogonal to the current dipolar 
model. Surprisingly, when analyzing the 5-year WMAP 
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data, iGroeneboom fc Eriksenl {2009) found supportive 
evidence for this model at the 3.8cr significance level, 
when considering angular scales up to i < 400. Thus, 
assuming that the WMAP observations are free of un- 
known systematics, there appears to be increasing evi- 
dence for both dipolar and quadrupolar structure in the 
CMB power distribution, at all angular scales. 

In this paper, we repeat the Bayesian analysis of 
lEriksen et al.l (|2007bf ). but double the angular resolution 
of the data. Nevertheless, we are still limited to relatively 
low angular resolutions, since the method inherently re- 
lies on brute-force evaluation of a pixel-based likelihood, 
and therefore scales as 0(N^ ix ). Yet, simply by spending 
more computer resources we are able to increase the pixel 
resolution from -/V s ido = 16 to 32 and decrease the degra- 
dation smoothing scale from 9° to 4.5° FWHM. This pro- 
vides additional support for multipoles between £ w 40 
and 80. While not sufficient t o provide a ful l and direct 
comparison with the results of lHansen et alj (|2008t ) , it is 
a significant improve ment over the results presented by 
lEriksen et all (|2007b[ ). 

2. OVERVIEW OF MODEL AND ALGORITHMS 

The Bayesian analysis framewo rk used in this paper ar e 
very similar to that employed by lEriksen et all (|2007bf ) . 
We therefore only give a brief overview of its main fea- 
tures here, and refer the reader interested in the full de- 
tails to the original paper and references therein. 

2.1. Data model and likelihood 

The starting point for our analysis is th e phenomeno- 
logica l CMB signal model first proposed by I Gordon et al.l 

d(n) = [l + /(n)]s(n) +n(n). (1) 

Here d(n) denotes the observed data in direction n, s(n) 
is an intrinsically isotropic and Gaussian random field 
with power spectrum Ce, f(n) is an auxiliary modulating 
field, and u(n) denotes instrumental noise. 

Obviously, if / = 0, one recovers the standard isotropic 
model. However, we are interested in a possible hemi- 
spherical asymmetry, and we therefore parametrize the 
modulation field in terms of a dipole with a free ampli- 
tude A and a preferred direction p, 



f(h) = A(h-p). 



(2) 



The modulated signal component is thus an anisotropic, 
but still Gaussian, random field, with covariance matrix 

S mod (n,m) = [l+A(n-p)]S iao (h,m)[l+A(m-p)], (3) 

where 

S iso (n» = rE( M+1 )W*" 1 )- 

47T 



(4) 



We now introduce one new fe ature compared to the 
analysis of lEriksen et al.l (|200 7b), for two reasons. First, 
we are interested in studying the behaviour of the mod- 
ulation field as a function of Grange, and therefore want 
a mechanism to restrict the impact of the modulation 
parameters in harmonic space. Second, we also want to 
minimize the impact of the arbitrary regularization noise 
(see Section[3]) on the modulation parameters at high ts. 
Therefore, we split the signal covariance matrix into two 



parts, one modulated low-£ part and one isotropic high-i? 
part, 

^total ^raod + S iso , (5) 

where only multipoles between 2 < I < £ mo d are included 
in S mo d, and only multipoles at £ > £ mo d are included in 
S; so . (Note that we are not proposing a physical mech- 
anism for generating the modulation field in this paper, 
but only attempt to characterize its properties. This split 
may or may not be physically well-motivated, but it does 
serve a useful purpose in the present paper as it allows 
us to study the scale dependence of the modulation field 
in a controlled manner.) 

Including instrumental noise and possible foreground 
contamination, the full data covariance matrix reads 



C = S nwd (A,p) + S is 



N 



(6) 



The noise and foreground covariance matrices depend on 
the data processing, and will be described in greater de- 
tail in g3 

We also have to parametrize the power spectrum 
for the under l ying is otropic component, Ci. Following 
lEriksen et all f2007bf) . we choose a simple two-parameter 
model with a free amplitude q and tilt n for this purpose, 



C t = q 



a 



lid 



(7) 



Here £q is a pivot multipole and Cf d is a fiducial model, 
in the following chosen to be the best -fit ACDM power 
law spectrum of iKomatsu et all (|2009f ). 

Since both the signal and noise are assumed to be 
Gaussian, the log-likelihood now reads 

- 2 log C(A, p, q, n) = d T C M + log | C| , (8) 

up to an irrelevant constant, with C = C(A,p, q, n). 

2.2. The posterior distribution and Bayesian evidence 

The posterior distribution for our model is given by 
Bayes' theorem, 

P(q,n,A,p\d,H) = 



C(q,n,A,p)P(q,n,A,p\H) 
P(d\H) 



(9) 



Here P(q,n, A,p\H) is a prior, and P(d|H) is a normal- 
ization factor often called the "Bayesian evidence" . Note 
that we now have included an explicit reference to the 
hypothesis (or model), H, in all factors, as we will in 
the following compare two different hypotheses, namely 
"HI: The universe is isotropic (A — 0)" versus "H2: The 
universe is anisotropic (A ^ 0)" . 

We adopt uniform priors for all priors in the follow- 
ing. Specifically, we adoptP(g) = Uniform [0.5, 1.5] and 
P(n) = Uniform[— 0.5, 0.5] for the power spectrum, and 
a uniform prior over the sphere for the preferred axis, p. 
The modulation amplitude prior is chosen uniform over 
[0, A max ], where A max = 0.15 is sufficiently large to fully 
encompass the non-zero parts of the likelihood. If more 
liberal priors are desired, the interested reader can eas- 
ily calculate the corresponding evidence by subtracting 
the logarithm of the volume expansion factor from the 
results quoted in this paper. 

With these definitions and priors, the posterior dis- 
tribution, P(q, n, A, p|d, H) is mapped out with a stan- 
dard MCMC sampler. The Bayesian evidence, E — 
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P(d\H ) is computed with the "nested sa mpling" algo- 
rithm (|Skillindl20dllMukherjee et alJl2006f >. For further 
de tails on both p r ocedur es, we refer the interested reader 
to lEriksen etall (|2007bT ). 

For easy reference, we recall Jeffreys' inte rpretational 
scale for the Bayesian evidence (Jeffr eys! 196 If ) : A value of 
A In E < 1 indicates a result "not worth more than a bare 
mentioning"; a value of 1 < A\nE < 2.5 is considered 
as "significant" evidence; a value of 2.5 < Ami? < 5 is 
considered "strong to very strong" ; and A In E > 5 ranks 
as "decisive". 

3. DATA 

In this paper we analyze several downgraded versions 
of the five-year WMAP temperature sky maps, namely 
the template-corrected Q-, V- and lU-band maps, as well 
as the "fore ground cleaned" I nternal Linear Combination 
(ILC) map (jGold et al.l l 2009h. Each map is down graded 
to low resolution as follows (|Eriksen et al.ll2007ah : First, 
each map is downgraded to a HEALPix 8 resolution of 
Aside = 32, by smoothing to an effective resolution of 4.5° 
FWHM and properly taking into account the respective 
pixel windows. We then add uniform Gaussian noise of 
<r„ = 1 [iK RMS to each pixel, in order to regularize 
the pixel-pixel covariance matrix at small angular scales. 
The resulting maps have a signal-to-noise ratio of unity 
at £ = 80, and are strongly noise dominated at ^ max = 95. 

Two different sky cuts are used in the analy ses, both of 
which are based on the WMAP KQ85 mask ([Gold et all 
2009). In the first case, we directly downgrade the KQ85 
cut to the appropriate A s jd e , by excluding any HEALPix 
pixel for which more than half of the corresponding sub- 
pixels are missing. This mask is simply denoted KQ85. 
In the second case, we smooth the mask image (consisting 
of 0's and l's) with a beam of 4.5° FWHM, and reject all 
pixels with a value less than 0.99. We call this expanded 
mask KQ85e. The two mask remove 16.3% and 26.9% of 
the pixels, respectively. 

The instrumental signal-to-noise ratio of the WMAP 
data is very high at large angular scales, at about 150 for 
the V-band at £ = 100. The only important noise contri- 
bution in the downgraded sky maps is therefore the uni- 
form regularization noise, which is not subject to the ad- 
ditional beam smoothing. We therefore approximate the 



noise covariance matrix by 



Note that this 



approx imation was explicitly validated bv lEriksen et all 
(|2007rD for the 3-year WMAP data, which have higher 
instrumental noise than the 5-year data. 

We also marginalize over a fixed set of "foreground 
templates" , t j , by adding an additional term to the data 
covariance matrix on the form Fj = ctititj, with > 
10 3 , for each template. In addition to one monopole 
and three dipole templates 9 , we use the V-ILC difference 
map as a template for both the U-band and ILC maps, 
the Q-ILC difference for the Q-band, and the W-ILC 
difference for the W^-band. However, these foreground 
templates do not affect the results noticeably in either 
case, due to the sky cuts used. 

4. RESULTS 



8 http://healpix.jpl.nasa.gov 

9 For an explicit demonstration of the importance of monopole 
and dipole marginalization on this specific problem, see Gordon 
2007a 
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Fig. 1. — Posterior distributions for the dipole modulation am- 
plitude, marginalized over direction and CMB power spectrum, 
computed for the KQ85 sky cut and £ moc i = 64. 



The main results from the analysis outlined above are 
summarized in Table [T] We consider nine different data 
combinations (ie., frequency bands, masks and multipole 
range), and show 1) the best-fit modulation axis and am- 
plitude, both with 68% confidence regions; 2) the sta- 
tistical significance of the corresponding amplitude (ie., 
A/ a a)', and 3) the raw improvement in x 2 and Bayesian 
log-evidence for the modulated model over the isotropic 
model. The last items are shown for the ILC with the 
KQ85 sky cut only. For r eference, we also quot e the ILC 
res ult for the Kp2 mask ((Bennett et al.ll2003bl ) reported 
by lEriksen et al.l (|2007b( ) when analysing the A s jd c = 16 
and 9° FWHM data. 

The reason for providing the full evidence for only one 
data set is solely computational. The total CPU cost for 
the full set of computations presented here was ~ 50 000 
CPU hours, and the evidence calculation constitutes a 
significant fraction of this. On the other hand, the ev- 
idence is closely related to the significance level Aj a a, 
and one can therefore estimate the evidence level for 
other cases in Table [T] given the two explicit evidence 
values and significances. We have therefore chosen to 
spend our available CPU time on more MCMC posterior 
analyses, rather than on more evidence computations. 

We consider first the results for the ILC map with the 
KQ85 mask and £ mo d = 64. In this case, the best-fit 
amplitude is A — 0.073 ± 0.022, non-zero at the 3.3<r 
confidence level. The best-fit axis points towards Galac- 
tic coordinates (l,b) = (224°, -22°, with a 68% uncer- 
tainty of 24°. The se results are consist ent with the 
results presented by lEriksen et al.l (|2007bf ). who found 
an amplitude of A = 0.11 ± 0.04 and a best-fit axis of 
(l,b)= (225°, -27°) for £< 40. 

Second, we see that these results are only weakly de- 
pendent on frequency, as both the V-band and W-band 
for the same mask and ^-range have amplitudes within 
0.5cr of the ILC map, with A = 0.080 and A = 0.074 and 
non-zero at 3.8<7 and 3.5er, respectively. (We have not in- 
cluded the Q-band analysis for the KQ85 mask, as there 
were clearly visible foreground residuals outside the mask 
for this case.) The corresponding marginal posteriors are 
shown in Figure [TJ clearly demonstrating the consistency 
between data sets. Figure [2] compares the best-fit axes of 
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TABLE 1 

Summary statistics for modulated CMB model posteriors 



Data 



Mask 



Significance (cr) A log C A log E 



ILC 

V-band 

V-band 

V-band 

W-band 

ILC 

Q-band 

V-band 

W-band 

ILC a 



KQ85 

KQ85 

KQ85 

KQ85 

KQ85 

KQ85e 

KQ85c 

KQ85e 

KQ85c 

Kp2 



64 
64 
10 
80 
(51 
64 
64 
64 
(51 
~ 40 



(224°,- 
(232°,- 
(224° , - 
(235°,- 
(232°,- 
(215°,- 
(245°,- 
(228°,- 
(226°,- 
(225°,- 



-22°) 
-22°) 
-22°) 
-17°) 
-22°) 
-19°) 
-21°) 
-18°) 
-19°) 
-27°) 



± 24° 
±23° 
±24° 
±22° 
±24° 
±28° 
±23° 
±28° 
±31° 



0.072 
0.080 
0.119 
0.070 
0.074 
0.066 
0.088 
0.067 
0.061 
0.11 



± 0.022 
±0.021 
± 0.034 
±0.019 
±0.021 
± 0.025 
± 0.022 
± 0.025 
± 0.025 
± 0.04 



3.3 
3.8 
3.5 
3.7 
3.5 
2.6 
3.9 
2.7 
2.5 
2.8 




Note. — Listed quantities are (first column) data set (first column); mask (second column); 
maximum multipole used for modulation covariance matrix, £ mo d (third column; marginal best- 
fit dipole axis (fourth column) and amplitude (fifth column) with 68% confidence regions indi- 
cated; statistical significance of non-zero detection of A (sixth column); the change in maximum 
likelihood between modulated and isotropic models, A log C = log £ mo d — log £i so (seventh col- 
umn); and the Bayesian evidence difference, A log E = log i? mo d — log Ei so (eighth column). The 
latter two were only computed for one data set, due to a high computational cost. However, 
other values can be estimated by comparing the significances indicated in the sixth column. 
a Results computed from A^ide = 16 and 9° FWHM data, as presented bv lEriksen et al.l (|2007bD . 




Fig. 2. — Posterior distribution for the dipole modulation axis, 
shown for the V-band map and KQ85 sky cut, marginalized over 
power spectrum and amplitude parameters. Grey sky pixels in- 
dicate pixels outside t he 2a confiden c e regio n. The dots indicate 
the axis 1) reported by Erikscn ct al. (2004a) in white; 2) for both 
the ILC and V-band maps (these have the same best-fit axis) with 
the KQ85 sky cut in blac k; 3) f or the tV-bands in blue, and the 
axis reported by Erikscn ct al. (2007b) in green. Note that the 
background distribution has been smoothed for plotting purposes 
to reduce visual Monte Carlo noise. 



th e three data se t s, and also indicates the axes rep orted 
bv lEriksen et alj (|2004aD and lEriksen et all (|2007bD . 

Next, we also see that the results are not strongly de- 
pendent on the choice of mask, as the amplitudes for the 
extended KQ85e mask are consistent with the KQ85 re- 
sults, even though it removes an additional 10% of the 
sky. However, we do see, as expected, that the error 
bars increase somewhat by removing the additional part 
of the sky, and this reduces the absolute significances 
somewhat. 

Finally, the best-fit modulation amplitudes for the V- 
band data and KQ85 mask are A = 0.11 for £ mo d = 40, 
A = 0.075 for £ mod = 64 and A = 0.066 for £ mod = 80 at 
3.5<7, 3.8cr and 3.7<7, respectively. This is an interesting 
observation for theoreticians who are interested in con- 
structing a fundamental model for the effect: Taken at 
face value, these amplitudes could indica te a non-scale 
invari ant behaviour of A 7 as also noted bv lHansen et al] 
(2008). On the other hand, the statistical significance of 
this statement is so far quite low, as a single common 



value A ~ 0.07 is also consistent with all measurements. 
Better measurements at higher €s are required to unam- 
biguously settle this question. 

5. CONCLUSIONS 

Shortly follo wing the release of t he first-year WMAP 
data in 2003, lEriksen etaH (|2004bD presented the early 
evidence for a dipolar distribution of power in the 
CMB temperature anisotropy sky, considering only 
the large angular scale s of th e WMAP data. Next, 
iGroeneboom &; Eriksenl (|2009f ) presented the evidence 
for a quadrupolar distribution of CMB power, and found 
that this fea ture extended over a ll £'s under considera- 
tion. Finally, Hansen et al. (2008) found that the dipolar 
CMB power distribution is present also at high £'s. The 
evidence for violation of statistical isotropy in the CMB 
field is currently increasing rapidly, and the significance 
of these detections are approaching 4a. 

In this paper, we revisit the high-£ claims of 
I Hansen et al.1 (|2008f) . by applying an optimal Bayesian 
framework based on a parametric modulated CMB model 
to the WMAP data at higher multipoles than previ- 
ously considered wi th this method, alb eit lower than 
those considered by lHansen et al.l ((2008). In doing so, 
we find result s very consistent with those presented by 
lHansen et al.l (|2008l ): The evidence for a dipolar dis- 
tribution of power in the WMAP data increases with 
£. For example, when considering the V-band data and 
KQ85 sky cut, the statistical significance of the modu- 
lated model increases from 3.2cr at £ mo d = 40, to 3.8<7 at 
^mod = 64, and 3.7cr at £ mo d — 80. 

The Bayesian evidence now also ranking within the 
"strong to very strong" category on Jeffreys' scale. How- 
ever, it should be noted that the Bayesian evidence is by 
nature strongly prior dependent, and if we had chosen a 
prior twice as large as the one actually used, the corre- 
sponding log-evidence for the ILC map would have fallen 
from Ami? = 2.6 to 1.7, ranking only as "substantial" 
evidence. For this reason, it is in many respects easier 
to attach a firm statistical interpretation to the posterior 
distribution than the Bayesian evidence. 

It is interesting to note that the absolute amplitude A 
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may show hints of decreasing with L It is premature to 
say whether this is due simply to a statistical fluctuation, 
or whether it might point toward a non-scale invariant 
underlying physical effect, in which case the amplitude 
A should be replaced with a function A{1). Either case 
is currently allowed by the data. 

To answer this question, and further constrain the 
overall model, better algorithms are required. The 
current approach relies on brute-force inversion of an 
iVpi x x iVpjx covariance matrix, and therefore scales as 
0(iVp ix ) or 0(N® ide ). However, already the present anal- 
ysis, performed at iV s ;d c = 32, required ^50 000 CPU 
hours, and increasing iV s i<ie by an additional factor of 
two would require ~3 million CPU hours. More efficient 
algorithms are clearly needed. 

To summarize, there is currently substantial evidence 
for both dipolar (Hansen et al. 2008 and this work) and 
quad rupolar power distribution (|Groeneboom fe Eriksenl 
|2009() in the WMAP data, and this is seen at all probed 
scales. The magnitude of the dipolar mode is consider- 
ably stronger than the quadrupolar mode, as a ~ 3.5a 
significance level is reached already at I ~ 64 for the 
dipole, while the same significance was obtained at £ ~ 
400 for the quadrupole. 

These observations may prove useful for theo- 
rists attempting to construct alternative mod- 
els for these features, either phenomenological 
or fundamental. Considerable eff orts have gone 
towards this goal already (e.g., lAckerman et al.1 



20071 iBohmer fc Motal [20081: [Carroll et all l2008allbl; 



Chang et alJl2008t lErickcek et al.ll2008allbl; iGordon et al 



2008allbl: 



Kahniashvili et al 



20051 _ iGumriikcuoglu et alj l2007t iHimmetoglu etal 



2008; 



Kanno et al 



20081: iKoivisto fc Motal l2008afbl iPereira et all 120071: 



Pitrou et al.1 l2008t iPullen fc KamionkowskH 120071 : 



Watanabe et al. I 120091: [Yokoya ma fc Sodal I2008T) . but 



so far no fully convincing model has been established. 
Clearly, more work is needed on both the theoretical and 
observational side of this issue. Fortunately, it is now 
only a few years until Planck will open up a whole new 
window on these issues by producing high-sensitivity 
maps of the CMB polarization, as well as measuring the 
temperature fluctuations to arc- minute scales. We will 
then be able to measure the properties of the dipole, 
quadrupole and, possibly, higher-order modes of the 
modulation field to unprecedented accuracy. 
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